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Summary 

Incidence of wave breaking for pure wind driven waves has been studied on Lake Washington at 
wind speeds up to 8 m s~ l . Video recordings were employed to identify and categorize the breaking 
events in terms ot micro-scale, spilling and plunging breakers. These events were correlated with 
the magnitude of the wave spectrum measured with a resistance wire wave gauge and band pass fil- 
tered between 6 and i 0 Hi. An equivalent percentage of breaking crests were found for spillinir and 
plunging events. v 6 

Wave forcing as measured by wind stress (or friction velocity, squared) and by inverse wave age, 

iu/C-, where C- is the phase velocity of the waves at the peak of the frequency spectrum, were 
found to be good predictors of percentage of breaking crests. When combined in a two parameter 
regression, those nvo variables gave small standard deviation and had a high correlation coefficient 
(66%). The combination of uz and u*/C p can be understood in physical terms. Furthermore, for the 
larger values of the dependence of wave breaking on wave age was stronger than at the low end of 

the values for u „ and u ./ C p , Thus, both the level of wave development as determined by inverse 
wave age, which we may term relative wind effectiveness for wave forcing and the wind forcing on 
the water surface determine the incidence of wave breaking. 

Substituting (J -j Q 75 (which is the dependence of whitecap cover found by Monahan and coworkere) 
anequivalent correlation was found to the prediction by it;. Slightly better standard deviation value 
and higher correlation coefficient were found by using a Reynolds number as predictor. A two-pa- 
rameter regression involving u; and a Reynold's number proposed by Toba and his colleagues (eg., 
Toba and Koga, 1986) which relates uz and peak wave frequency, improves the correlation even 
more but is less easy to interpret in physical terms. 

The equivalent percentage of breaking crests obtained in our previous study (Weissmanet al., 1984) 
wasreported at 8.6% for a short record obtained at U l0N of about 6 ms~^. Typical values in the airrent 
study for similar conditions are 6%, which is consistent with the previous study in view of the scatter 
In that study we did not have a video recording system, so the observed breaking may include more of 
the micro-scale breaking events, and the value, 8.6%, is well within the range of highly probable 


INTRODUCTION 

Wave breaking can be t dramatic change in sea surface characteristics and occurs on many scales. It 

plays a significant role in air- sea exchanges of momentum, energy and mass in various ways (e.g., 
Donelan, 1990; Banner, 1990). In particular, it is the dominant mechanism responsible for wave dis- 
sipation; it removes momentum and energy from the wave field and transfers the momentum to sur- 
face currents and the energy into both turbulence and currents. The dissipation of wave energy, due 
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to breaking, is one of the ke y parameters in wave prediction models which are based on the energy 
transfer equation: 

i£ + v-VE-s in + s„, + s*, (i) 

where £ is the wave spectral energy density, V is velocity at which the wave energy propagates (i.e., 
group velocity plus currents), S tn is source function due to wind input, 5„/ is source function due to 
nonlinear wave-wave interactions, and Sdiu is source term due to dissipation. In the state of the art, 
third generation ocean wave prediction model, the WAM model (the WAMDI Group, 1988), S tn 
based on field measurements has been adopted from that of Snyder et ai. ( 1 98 1 ), derived on theo- 

retical grounds has been taken from Hasselmann et al. (1985) and, inferred from the residual 
term in numerical energy balance experiments has the form proposed by Komen et ai. (1984). De- 
spite its importance. is the least known of these three source terms. 

Whitecaps, which are related to the wave breaking but also includes the more persistent visible foam 
patch, have been studied extensively by Monahan and coworkers using statistical analysis of photo- 
graphs and video recordings (e.g., Monahan, 1968, 1971; Monahan and O’Muircheaxtaigh, 1980; 
Monahan et al., 1983; O’Muircheanaigh and Monahan, 1986; O’Muircheartaigh et al., 1991). In 
these studies areal whitecap coverage has been found to correlate with the 10 m neutral wind speed 
raised to the power 3.75, (also sec, Wu, 1979, 1986). 

Because of the increased significance of small scale perturbations on the sea surface for the return of 
incident radar signals and for microwave emissions, renewed interest in defining wave breaking and 
relating its occurrence to other measurable quantities has developed (e.g.. Banner and Fooks, 1985). 
Jessup etaL (1990, 1991a) and Bushel al. (1991) have related the sea spikes in radar return at X,C 
and Ku bands to wave breaking observed from video records. The latter studies found that sea spikes 
often occurred without visible breaking events, but Jessup et ai. (1991b) comment that the largest 
spikes in their study were all directly relatable to readily visible plunging breakers (i.e. , dramatically 
breaking wave crests). 

The Bragg scatter is the dominant mechanism responsible for the variations in the normalized radar 
backscaner cross section, o 0% as a function of wind speed or stress. It is currently recognized that 
breaking events may modify this relationship, especially since it is unlikely that occurrence of wave 
breaking has the same dependence on wind speed or stress as the part of the wave spectrum corre- 
sponding to Bragg scattering water waves. It has also been observed that wave breaking is related 
not only to the wind forcing but that details of the underlying sea state must also be considered. Nu- 
merous studies have considered the hydrodynamic instability of the water surface itself, i.e., in situa- 
tions with no atmospheric forcing at all (see other articles in this volume). On theoretical grounds 
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Phillips (1988) predicted that contributions to o 0 by breaking waves should be proportional to uj, 
where h* is friction velocity defined through r 0 /q * u\. x 0 is wind stress atthe sea surface and Q is 
air density. Jessup et ai. (1990, 1991b) verified this relationship. 

Thus, the variability of processes responsible for wave breaking events may also be responsible for 
some of the scatter found between measured radar cross sections and wind speed (or wind stress ). So 
far as the relationship to wind stress is concerned, it is quite possible that the intensity of wave break- 
ing affects the wind stress in a manner which more closely parallels the effect on radar backsc alter 
than the relationship between 10 m neutral wind speed, U and <7 0 . It is customary to use U 
rather than the tree wind speed in this relationship. U IW is convertible to a wind stress by defining a 
relationship between the neutral drag coefficient, C^ f and U viz: 

T * “ Q Cdn U\qn • ® 

This relationship represents an average condition if typical values for C DN are used (e.g.. Smith, 
1980; Liu et al., 1979). For anomalous wave breaking a different C DN value should probably be 
employed and <7* may be a very good measure of the effective r 0 . 

In this study we extend the work of Weissman et al. (1984), who developed a technique for identify- 
ing breaking waves from the record of a resistance wire wave gauge. That study was based on a lim- 
ited time series and only visual observations of the conditions at the wave wire. In this study video 
recordings were employed. 

A preliminary study on sea spikes and wave breaking by Bush et al. (1991 ) with data obtained on our 
fak e Washington site has been carried out. Further work on the relationship between wind and wave 
parameters, frequency of occurrence of wave breaking and radar cross-sections at X, C and Ku 
bands will be repotted elsewhere in the near future. 

EXPERIMENTAL SET-UP AND DATA SETS ANALYZED 

This work was earned out at our field station on Lake Washington where wind and waves arrive after 
an over water fetch of 7 km (Figure 1 .a). A complete suite of environmental measurements was col- 
lected, which includes mean wind speed, temperature and humidity at 2 m and 8 m heights, turbulent 
fluxes of momentum, heat and water vapor at 8 m and wave height from a resistance wire wave gauge 
hanging from a boom extending 2 m upwind of the supporting platform (Figure 1 .b). The wave wire 
was made of stainless steel, 100 \im in diameter. 

For the purpose of this study, 10 data sea were analyzed. The length of each data set is 2 hoots. The 
experimental conditions during these data sea are summarized in Table I where for each run the 
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Figure la: The location of our field station 
(MS MAST) on Lake Washington. The tower 
is 15 m offshore and at a depth of 4 m. The 
contours show the water depth in meters. 
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Figure lb: The Lake Washington tower with 
meteorological instrumentation, wire wave 
gauge (at the left most tip of the low level 
boom) and the video camera (at the lower 
section of the platform). 


mean values (the first line) and the standard deviations (the second line) of various environmental 
variables are provided. 

ANALYSIS 

Weissman et ah ( 1984) defined three types of breaking events; micro-scale, spilling and plunging. 
The three types are fairly easily separated. Typical examples are seen in Figure 2. Two independent 
observers of the video record agreed on the classification of the breaking events into these three cate- 
gories with insignificant deviations. For the event to be counted in the video record it had to occur 
within aO.l m radius of the penetration of the wave wire through the mean water level. 

For the wave breaking analysis, wave wire data were processed to obtain a time series of the spectral 
energy density, f l0 , in the frequency band 6-10 Hz with a sampling interval of 1/8 second (see 
Ataktiirfc, 1991). In our previous study (Weissman et ah, 1984) we had devised a technique to detect 
the breaking events from such a time senes by using the breaking criteria: (i) the spectral energy in a 






123 


high frequency band exceeds a threshold, and (ii) the data point is in the crest region (i.e.. in the 
vicinity of a local maximum). Since the threshold, i.e., the absolute value of the measured band 
energy vanes with wind speed (stress), gustiness, underlying long waves, currents etc., it had to be 


determined individually for each run. In the current study, this difficulty was greatly reduced by de- 
fining a new detection parameter, N„, based on the fluctuating component of band energy normal- 
ized by its root— mean— square value; 


N, 


^S-IO Z«- 


10 


"<7,6- 10 


(3) 


Tabic I: General description of the data sets analyzed. Duration of each run is two 
h°un. For each run, the first row indicates the mean values and the second row the 
standard deviations. 
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Y/M/D 

Time 
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<°C) 
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z/L 

fp 

{Hi) 

H, 

(m) 

ak 


1100 
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SOI 

0.16 

-1.15 

2336 

-0.27 

0.67 

IBM 

0.09 

1.88 j 




ESI 

0.02 

0.22 

0.03 

0.04 
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IKB 

0.01 

0.23 
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86/8/25 
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(a) micro— scale (b) spilling (c) plunging 


Figure 2: The three types of wave breaking characterized from our video records. 


where £ 6 _ l0 and 10 respectively, arc the mean and standard deviation of £ 6 _ l0 calculated 
from a record typically 1 7 minutes in length. (Out of 66 records, two records were short -8.5 minute 
in duration. The statistics from these runs were included by simply doubling the number of breaking 
events.) Also, data points within ±90* phase of the local maximum were considered to be on the crest . 

Long wave spectra were calculated from the 17 min wave height records. An equivalent number of 
wave crests, N crt „ , in the record was defined as the duration, T , times the frequency of the waves at 
the peak of the spectrum. f p: 


Her*,, « Tf p , 


(4) 


From video records, N br€0kt number of crests on which breaking occurred, was found. Percentage of 
crests with breaking, 




N, 


hrro* 




x 100 


(5) 


was initially calculated separately for the events of plunging, spilling and micro-scale breakers. 
Since the frequencies of occurrence of spilling and plunging breakers were about the same and the 
micro-scale breaking turned out to be ubiquitous and not well correlated with any turbulence or sea 
state measure, presentation of results has been done in terms of the sum of percentages of crests with 
breaking of type plungers and spiilers, 


In a separate study, AtaktUik (1991) has analyzed the wind stress measurements of the current study 
and additional data runs, and obtained the relationship between the drag coefficient and neutral 10 m 
mean wind speed valid for the fetch and water surface conditions at our 1 jJce Washington site; 


125 


C DN = (0.75 +0. IQU IW ) x 10" 3 (6) 

This formula was used with the Liu et ai. (1979) formuiation for calculating turbulent fluxes from 
mean atm ospheric measurements. The calculation of wind stress or equivalently the friction velocity 
from mean wind speed with this numerical scheme includes corrections for the influence of atmo- 
spheric stratification. 

We also determined the inverse wave age, f/ lcw /Cp or u */Cp t for each run t where = g/ajpisthe 

phase speed of the dominant waves corresponding to the angular frequency, c o p = 2n of the spec- 
tral peak and, g is the acceleration due to gravity. 

cf 0 Q was correlated with various measures of the atmospheric turbulence and the relative wind 
forcing (identically equivalent to the inverse wave age). Both single and dual parameter linear re- 
gressions (Lapin, 1983) were calculated for several variables using the statistical program package 
Quattro-Pro. A curvilinear fit was also tried when it was noted that a linear fit for a certain variable 
was not optimal. 

RESULTS 

Figure 3a is a plot of %B determined from video records and separated into micro-scale, spilling and 
plunging types, versus (J . Micro-scale breaking occurs to about the same extent at all wind 
speeds, while spilling and plunging type breakers display a similar behaviour in frequency of occur- 
rence which increases with increasing U 10W . Therefore, further analyses were carried only for the 
combined effects of spilling and plunging breakers. 

Comparisons of simultaneous video and wave wire records showed that, in general, micro-scale 
breaking is associated with events for which N 0 < 8 and, spilling or plunging breaking may be char- 
acterized by N a £ 8. Events with N a a 12 were always due to plungers. In Figure 3b, percentage 
of crests with spilling or plunging breaking, %B s +p, determined from video records and percentage 
of crests with events N 0 ^ 8 obtained from wave wire records are plotted as a function oi U {W . 
Although the distributions of the data points are similar, some differences in their magnitudes are 
visible: values from video records are smaller at low winds, and larger at high winds. These differ- 
ences may result from, respectively, (i) sensitivity of wave wire to small scale features which may 
not be visible on video records, (ii) insensitivity of wave wire to events that occur nearby but do not 
directly touch the sensor. The latter may be improved by using an array of wave wires which would 
allow spatial averaging. 

In the studies by Jessup etal. (1990, 1991a) andBushetaL(1991)some spikes in the radar backscat- 
ter were found not to correspond to any visible breaking events. Those events may be carried by the 
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Wind speed, Ujon {mis} 


Figure 3a: Percenuges of crests with micro-scale, %B mt spilling, %B tf or plunging, %B p , 
breaking versus wind speed, Ujon, determined from video recor ds . 


same features (for instance sharp comers) that cause N 9 to be large (but still < 12). This indi- 
cates that the radars are quite sensitive to such features, even if they do not cover the whole radar 
footprint. 

Some of our 8 < N c < 12 events (or perhaps even a few of the No > 12 cases) correspond to 
strong Doppler shifts causing the amplitudes of the signals at 6-10 Hz . to be increased (being that 
they are really due to surface features of much lower intrinsic frequency. See, for instance, Ataktiirk 
and Katsaros (1987) for a discussion of Doppler shifts in wave wire data.). 

Figure 4 shows the dependence of %B 9 +^ on/te, £/*J 3 , it*/C,and ui Reynolds number, was 
defined as (Toba and Koga, 1986): 

m j± 

“ VCJp 


Re 


(7) 
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Figure 3b: Statistics of breaking events obtained from video and wave wire records versus 
wind speed, Uiqn. 


where v is the kinematic viscosity of air. All these parameters are about equally good predictors of 

However, from the statistics of the regression analysis provided in Table Q it is seen that u\ 

which explains 62% of the variations in %£ I+r is the best one among the individual variables used. 
(Although ul* gives a slightly better correlation, for the limited data set the difference is not signifi- 
cant.) 

This point is further illustrated in Figure 3 of %B,+ p versus it./C, where the data points have been 
stratified according to the values of There is obviously a dependence on both the value of u \ and 

u ./ C p . Even though the two measures are not independent (since both contain it •), it makes sense to 
examine for a certain wave age the variability in the %B,+ p as a function of wind stress (here calcu- 
lated from U lQN ). The regression lines show that for strong relative wind forcing (large inverse wave 
age), the increase in %B,+ p as a function of u./C f is greater for large values of uj. Also, from Table 
II it is seen that the pair, u\ and u*/ C ^ explains a higher proportion of the variance of %B s + p than ui 
alone with a proportional reduction of 1 1 % in the previously unexplained variation. 

Multiple regressions using (ul.Re) or (u\ t Re % u./C p ) produced still better statistical results by ex- 
plaining, respectively, 69% and 72% of the variance of %B 1+ ^. However, physical interpretations of 
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Figure 4: 1 incar regressions of %B f +p for various predictors. A summary of the statistics ob- 
tained from regression analysis is provided in Table II. 


Table Q: Summary of the statistics obtained from regression analysis. The analyses 
were earned out following the definitions in Lapin (1983). Square-root of the 3rd 
column is the correlation coefficient. 


Regression 

Standard Error 
of Estimate 

Sample Coefficient 
of Determination 

=> 0.03 + 2.68 x I0" 3 rte 

1.07 

0.54 

0.44 + 2.32 X 10 _J l/$2J 

1.02 

0.38 

- - 3 . 97 + 79 . 48 u./C, 

0.99 

0.60 

- - 0 . 41 + 61 . 09 ui 

0.99 

0.62 

- - 2 . 47 + 32 . 87 «i + 42 . 37 ii./C, 

0.93 

0.66 


such combinations of v enables are not as clear to us (but see Toba and Koga (1986)). Also, in view of 
the large scatter, it does not make sense to try very complicated functional forms until a greater data 
base has been formed or a convincing physical relationship has been proposed. Therefore, these re- 
sults are not included in Table 0. 
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a. 04 0.06 0.06 0.10 


«vc. 

Figure 5: Percentage of crests with spilling or plunging breakers, %B s + pt versus «• /C p . 
For three ranges of u \ separate symbols and linear fits have been used, as indicated on the 
figure. 


DISCUSSIONS 

Predictability of the percentage of crests with spilling or plunging breakers, were investi- 

gated using a set of variables. Re, u m )C p and These variables were chosen because such 
correlations make sense from a physical point of view. (For interpretations of the first two variables, 
see Toba and Koga ( 1 986) and Wu (1979), respectively). In Table II we note that the highest sample 
coefficient of determination and the minimum standard error ate obtained for the regression of 
%fl J+F as a function of u, and u./Cp. This result implies that for wave breaking both the magnitude 
of the wind stress and the degree of evolution of the surface waves are important. 

It is easy to understand that young waves which are being strongly forced by a large wind stress 
should be breaking more frequently. Whether young waves* which are also short (i.e., the wave- 
length at the peak of the spectrum is short) break more firequentiy than those that have a wider spec- 
tral range through which to distribute the wind input by non-linear interactions - we intend to inves- 



tigtte further. We hypothesize that the absolute width of the wave spectmm (or the value of the wave 
number at the peak of the spectmm) could be considered as another parameter for improving the 
prediction of the incidence of wave breaking. 


Jessup et al. (1990) suggested that the information content of the frequency of spikes in the radar 
cross section may also be used to study the variability of wave breaking. We suggest that the band 
passed data from an array of simple wave gauge wires may serve the same purpose. Since we have 
seen that N 0 > 12 consistently corresponds to visible breaking in all our wind wave records, this 
characteristic can now be used to count breaking waves. Enough cases can then be added by objec- 
tive analysis to establish statistically significant information about the dissipation of surface waves 
under various forcing conditions as determined by wind stress and relative wind stress forcing. We 
are currently adding more cases on the breaking of pure wind waves on Lake Washington from two 
mote summer seasons. 


Once the interrelationships have been established for the simple case of wind driven waves, the role 
and importance of swell (of variable amplitude and relative direction to the wind waves) as well as 
surface currents in wave breaking can then be studied more effectively. 
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